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UTILISATION PRATIQUE 


1 Introduction 


L'objet d'HYDROLAB n'est certainement pas de remplacer des logiciels beaucoup plus complets tels 
qu'on peut en trouver dans le commerce. 

Notre volonté a été tout d'abord de créer un outil très simple parfaitement intégré au logiciel EXCEL 
dont l'utilisation est universelle. Pour se servir correctement d'HYDROLAB, il est indispensable d'avoir des 
bases sur l'utilisation d'EXCEL. 

Il est préférable d'utiliser HYDROLAB avec un environnement utilisant la virgule comme séparateur 
décimal. Cependant, on peut également utiliser HYDROLAB avec le point comme séparateur ; dans ce cas de 
figure il faut par contre remplacer manuellement les virgules par des points dans les sorties des programmes de 
régression multiple, d'analyse en composantes principales et de comblement des lacunes. Il suffit de sélectionner 
l'ensemble de la feuille puis d'utiliser la fonction "Remplacer" du menu "Edition". 

D'une façon générale, HYDROLAB est donc interfacé via EXCEL à la plus part des autres logiciels tels 
que : WORD, SURFER, ... 

Le second objectif est de répondre aux questions les plus fréquemment posées aux Hydrologues. Ces 
questions portent essentiellement sur : 

- L'analyse univariée (A justements) ; 

- L'analyse multivariée (Régressions multiples) ; 

- L'analyse en composantes principales (ACP) ; | Ces points sont traités par des macros commandes. 

- Le comblement de lacunes dans des séries de 
données ; 

- La détection d'anomalies dans les séries de données ; 
- L'analyse spatiale ( variographie ). 


- Des fonctions statistiques classiques telles que F de 
Fisher-Snedecor, intégrale de Gauss... ; 
- Des fonctions liées à l'estimation de l'évapo- 
transpiration potentielle (Penmann, Durée du jour, 
radiation, Mc culloch, FAO...) ;|Ces points sont traités par des "fonctions 

- Une fonction pour passer de l'évapotranspiration | personnalisées" au sens d'EXCEL. 
potentielle à la réelle ; 
- Des fonctions de passage des coordonnées 
géographiques à différentes coordonnées Lambert. 


Vous avez ici la version 98.2. Elle est l'aboutissement actuel d'une lignée de commandes écrites depuis 
1990 dans les versions 2.2 à 97 d'EXCEL. Au fur et à mesure, HYDROLAB s'enrichit, se perfectionne et se 
corrige en fonction des utilisations qui en sont faites. Aussi n'hésitez pas à nous faire part des anomalies que 
vous auriez détectées, des améliorations que vous proposeriez et des développements que vous souhaiteriez, à 
l'adresse suivante : 

UPRESA 6046 du C.N.RSS. 
98 Bd E. Herriot -B.P. 209 
06204 NICE Cedex 3 
Fax; 04 93 37 54 30 
E. mail :[laborde@unice fr 

Les données utilisées à titre d'exemple dans ce guide sont fournies dans classeur sous la forme des trois 
fichiers EXCEL suivants : 
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2 Les macros commandes 


2.1 Ajustements 
HYDROLAB permet de réaliser des ajustements aux lois les plus fréquemment utilisées en hydrologie : 


Loi de Gauss (loi normale) ; Loi de Fuller ; 
Loi de Galton (loi log-normale) ; Loi de Weibull ; 
Loi Racine-normale ; Loi de Poisson. 


Loi de Gumbel ; 


2.1.1 Règles communes : 


2.1.1.1 Vos données sont préservées 


D'une façon systématique, les feuilles Excel où sont situées vos données de départ ne sont jamais 
modifiées. Les feuilles de commandes et de présentation des résultats utilisés par HYDROLAB sont verrouillées. 
Vous pouvez les modifier mais vous ne pouvez pas enregistrer les modifications par erreur. Au cas ou vous 
auriez modifier une feuille originale d'HYDROLAB , il suffit de la fermer puis de la rouvrir. 


Si vous voulez conserver vos résultats pour les faire figurer dans un rapport, il suffit de copier ce qui vous 
intéresse dans les feuilles Excel puis de les copier sous Word. 


Si vous voulez enregistrer les résultats sous forme de fichier Excel, il convient d'utiliser la procédure 
"Enregistrer sous” du menu "Fichier". Vous pouvez alors choisir librement les noms des feuilles Excel à 
enregistrer. Attention ! Il est préférable d'enregistrer en premier les feuilles de graphique (extension .xlc) puis 
seulement les feuilles de calcul (extension .xls). 


2.1.1.2 Sélection de la plage de données 


L'échantillon que l'on veut ajuster doit être nécessairement une suite de valeurs numériques (ou des 
formules numériques) situées sur une seule et unique colonne d'une feuille EXCEL. Si ce n'était pas le cas vous 
devez vous y ramener par les différentes fonctionnalités d'EXCEL. 


Il suffit tout d'abord de sélectionner la plage de cellules à traiter. Sur la figure 8 on s'apprête à traiter les 
17 pluies annuelles à Cannes de 1959 à 1976. Sur la figure 9 nous avons une autre sélection tout à fait valide 
concernant des pluies à Antibes de 1960 à 1974 à l'exclusion de l'année 1968 pour laquelle il y' a une lacune. 
Pour obtenir cette sélection multiple, il a suffit de sélectionner la plage 1127 à 358 mm puis de maintenir 
enfoncée la touche Ctrl et sélectionner la plage 768 à 1075 mm. Par contre la sélection illustrée par la figure 10 
n'est pas valide puisqu'on a pris en compte un terme alphanumérique correspondant à la lacune d'observation de 
1962 à Clans. 
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Fig.8 : Sélection continue valide  Fig.9 : Sélection discontinue valide Fig. 10 : Sélection non valide 


2.1.1.3 Lancement de l'ajustement 


Une fois les données sélectionnées, il suffit pour 
procéder à un ajustement, de choisir dans la barre de 
menu HYDROLAB, le type de loi désiré. 


Dans l'exemple de la figure 11 les pluies à 


Cannes de 1959 à 1976 seront ajustées à une loi 
Normale. 


2.1.1.4 


Les feuilles de résultats 
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Fig. 12 : Exemple de résultats graphiques 


On trouve également des informations 
succinctes sur la taille de l'échantillon traité, les 
paramètres ajustés et l'intervalle de confiance tracé. 


On peut toujours habiller ce graphique en 
utilisant les fonctionnalités d'EXCEL En cliquant sur le 
titre ou les axes on peut en modifier les textes, les 
polices ... Par "affichage", puis "barre d'outil" on peut 
accéder à toute forme de traits, de motifs ... La figure 
13 en est un exemple. Il suffit lorsque le graphique 
vous satisfait, de le copier (Ctrl C) puis de la coller par 
exemple dans votre rapport sous WORD. 
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Fig.11 : Lancement d'un ajustement 


Quelques instants après le lancement, vous 
verrez apparaître les feuilles de résultats. Elles sont 
toujours de même type et superposées. Celle qui vous 
est visible, contient la représentation graphique de 
l'ajustement et les paramètres essentiels. Comme 
l'indique la figure 12, les feuilles de graphiques portent 
l'extension ".xlc" (ici précédée de Norm puisqu'il s'agit 
d'une loi de Gauss). 


Les valeurs de la variable sont toujours portées 

en ordonnées et les fréquences au non dépassement en 
abscisses, selon les échelles adaptées à la loi. 
La loi théorique est matérialisée par une droite en trait 
fort. La plus part du temps vous verrez également deux 
courbes en trait fin limitant un intervalle de confiance à 
80% ( intervalle que vous pourrez modifier). 


Fig. 13 : Exemple d'habillage personnalisé 


Une deuxième feuille est située sous le graphique, vous pouvez y accéder soit par le menu "fenêtre" où 
apparaît le nom d'une feuille de calcul avec le suffixe ".xls", soit en diminuant la taille de la fenêtre du 
graphique. On peut bien sur avoir les deux fenêtres visibles en mosaïque. 


Les premières colonnes sont toujours les 
mêmes, les valeurs dans l'ordre de départ, les valeurs 
Des Te | on [ten | TE] Mie | vue Lies | ee rangées par ordre croissant, l'ordre de classement, la 

= | ej | R | = = | | = | fréquence expérimentale calculée par la formule : 


RES a a EU né las = 


| ml à = | as | s | a (x; )= an (i rang de classement, n taille de 
= à: | i n 
z = f : E : - l'échantillon) 


= m | à = : a% En y trouve également, les paramètres de 
i ] = W | 1m l'ajustement, la taille de l'échantillon et l'intervalle de 
5 Û confiance exprimé en pour-cent. On peut modifier la 
valeur de 80% par défaut en remplaçant simplement la 
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Fig. 14 : Exemple de résultats numériques 


Dans le bas de la feuille, on trouve un second tableau qui permet d'obtenir des valeurs quantiliques ou au 
contraire trouver la fréquence et la période de retour d'une valeur quelconque. Il suffit de mettre la valeur 
souhaitée à la place des points d'interrogation. Dans le cas où l'on désirerait un tableau comportant plusieurs 
estimations de ce type, il suffit de sélectionner les cellules comprises entre celle sous "Fréquence" et celle sous 
"Per. Ret." , puis de les recopier vers le bas sur autant de lignes que nécessaire. 
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Dans le tableau ci-contre, 
l'intervalle de confiance a été ramené 
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période de retour 5, 10 et 20 ans, 
ainsi que les périodes de retour de 
pluies de 400, 800 et 1200 mm. 
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Fig. 15 : Modification d'une feuille 
de résultats numériques 


2.1.2 La loi de Gauss 
La mise en œuvre de la loi de Gauss correspond intégralement au dispositif commun de lancement et 

de présentation des résultats. Cette loi est définie sur H o0,+00] et peut convenir à n'importe quel échantillon. 
Il nous parait seulement bon de préciser quelques points techniques. La fonction de répartition 
u? 


Ea X—X 
e 2 duavec u= 


dépends de deux paramètres qui sont la moyenne X et l'écart-type ©. 


F(x) = rl 


e Ces paramètres moyenne et écart type sont évalués selon les estimateurs sans biais suivant : 


(paramètres estimés par la méthode des moments) 
e L'intégrale de Gauss est approchée par les formules suivantes : 


Calcul de u pour F <= 0,5 : 


-u=Xx- Erreur avec x = Erreur ! 


(pour F > 0,5 : on fait F = 1 - F et u = - u) 


Calcul de F pour u > 0 : 
0,5 


F (u)=1- 


(1+0,196854 x! + 0,115194 x? + 0,000344 x° + 0,019527 x* + 
(pour u < O on fera u=-uetF=1-F) 


e pour l'intervalle de confiance à œ % sur le quantile xp les bornes sont évaluées par : 
XF + Erreur !6 (Informations Techniques du CTGREF, Cahier 31, N°2, 1978) 
(2 fois plus pour la borne supérieure et 2 fois moins pour la borne inférieure) 
tF : variable réduite de Gauss ayant la fréquence au non-dépassement F 


; ga : f l-a 
ta : Variable réduite de Gauss ayant la fréquence au non-dépassement 1 — ES 


2.1.3 La loi de Galton ( loi Log-normale) 


La mise en œuvre de la loi de Galton correspond intégralement au dispositif commun de lancement et 
de présentation des résultats. Il faut cependant préciser quelques points techniques : 
u u? 
1 ne log(x—-x,)—a 
-= fe 2 du avec a 2e) 2 
V27 ” b 


co 


La fonction de répartition est F(x) = ; ici nous considérerons 


que c'est la variable log(x-x,) qui suit une loi de Gauss à deux paramètres : 


- a: la moyenne des log(x-x,) ; 
- _b:l'écart-type des log(x-x,) . 


D loga; —X0) 
_ i=l 


n n—-] 


Dog” (x —Xo)—n a? 
i=l 


a 


Dans cette approche on voit bien que le paramètre de position x, n'est pas considéré comme un 
paramètre d'ajustement mais comme une caractéristique connue a priori (borne inférieure de l'intervalle de 
définition). Par défaut, x, est supposé nul, mais vous pouvez lui donner une autre valeur en l'entrant directement 


dans la cellule située à droite de " x, =". Il faut évidemment que x, soit strictement inférieur à la plus petite 
valeur de l'échantillon. 
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Fig. 16 : Modification du paramètre de position dans une loi de Galton 
(à gauche x, est pris à priori égal à 0, à droite on s'est imposé x, = 5 en le tapant directement dans la feuille) 


Pour le calcul des quantiles et de leurs intervalles de confiance, les calculs sont faits sur la variable 
transformée log(x-xo) supposée gaussienne, puis ramenés en variable naturelle par x = Xo + 10 #70, 


2.1.4 La loi Racine-normale 


La mise en œuvre de la loi de Racine-normale correspond intégralement au dispositif commun de 
lancement et de présentation des résultats. Comme pour la loi de Galton on se ramène à une loi de Gauss par un 
changement de variable : 

e ?duavec u= 


La fonction de répartition est F(x) = ; ici nous considérerons que 


2 
1 f — x —a 
V27 * b 
c'est la variable x qui suit une loi de Gauss à deux paramètres : 
- a: la moyenne des Jx: 


- b: l'écart-type des Jx à 


Cette loi est définie sur [+0] et ne peut donc convenir qu'à une population de valeurs positives. Pour 
le calcul des quantiles et de leurs intervalles de confiance, les calculs sont faits sur la variable transformée vx 


; ; : g . T3 
supposée gaussienne, puis ramenés en variable naturelle par x =4Xx ^. 


2.1.5 La loi de Gumbel 


La mise en œuvre de la loi de Gumbel correspond intégralement au dispositif commun de lancement et 
de présentation des résultats. La fonction de répartition est : 


X-X 
F(x)=e® avec u= g 


8 
Le terme u est la variable réduite de Gumbel ; xQ est le paramètre de position (mode) et g est le 


(variable réduite de Gumbel) 


paramètre d'échelle différent de zéro et positif (g est aussi appelé "gradex"). 
Cette loi est définie sur [- 00,00] et peut convenir à n'importe quel échantillon. 


On peut noter dès à présent le comportement asymptotiquement exponentiel de la distribution de Gumbel : si 
F tend vers 1, en posant T = Erreur !, la variable réduite de Gumbel u tend vers Ln T. 


Les paramètres sont estimés par la méthode des moments : 


g=0,786et x= x, -0,577 s 


L'intervalle de confiance à & % sur un quantile xp s'exprime en fonction de l'écart-type © par : 
ĉr -h] G<xF< ĉr +h 0 
où h] et h2 sont des paramètres dépendant de la taille n de l'échantillon de la fréquence F et de la valeur de ©. 
e hı et h? seront évalués par la formule suivante (avec le signe + pour h2 et le signe - pour h1) : 


2 
Ba J1+113tp +11t} +7% (Lite +0,57) 
n 


h1,2 = 2 3 (Informations Techniques du CTGREF, Cahier 31, N°2, 1978) 


1-117% 


n 


’ ie 3 5 ; l-a 
° u4 est la variable réduite de Gauss correspondant à la fréquence au non-dépassement 1- —— 


e tp est la variable réduite de Gumbel correspondant à la fréquence au non-dépassement F, ramenée à sa 
moyenne et à son écart-type : 


tp = Erreur ! 


2.1.6 La loi de Fuller 


La mise en œuvre de la loi de Fuller correspond intégralement au dispositif commun de lancement et de 
présentation des résultats. La fonction de répartition est : 


+ #0, 
F(xX)=1-e £ 
Cette loi a donc deux paramètres (x, et g) et est définie sur l'intervalle [xo;+o]. Nous considérons ici 
que x, est la borne inférieure de l'intervalle de définition et qu'il est donc connu a priori. Dans la feuille la 
valeur par défaut de x, est la plus petite valeur de l'échantillon mais on peut toujours se donner une autre valeur 


de x,en la tapant directement dans la cellule située à droite de " x;=". Il faut bien sûr, que x, soit inférieur ou égal 
à la plus petite valeur de l'échantillon. 
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Fig. 17 : Modification du paramètre de position dans une loi de Fuller 
à gauche x, est pris à priori égal à 4(plus petite valeur de l'échantillon), 
à droite on s'est imposé x, =0 en le tapant directement dans la feuille 


Le paramètre g (que l'on peut également appeler gradex) est estimé par : g = X — X9 
L'intervalle de confiance à & % sur un quantile xp s'exprime ainsi : 


2 ü (Rr —xo) (Rr —xo) 


"<< xp <Èp+u 
a Jn F F œ Jao 


; 2e Nr ; 1-4 
e u4 est la variable réduite de Gauss correspondant à la fréquence au non-dépassement 1- —— 


(formule simplifiée d'après MIQUEL 1984 et valable si g est le seul paramètre ajusté). 


2.1.7 La loi de Weibull 


La fonction de répartition de Weibull est: 
F (x) = 1 - eP X - Xo)P 
… dont la densité est : 
KQ =P p @-xo)P"l e-P © - xo)P 
Dans cette loi figurent trois paramètres p, p et Xọ, mais ce dernier, borne inférieure de l'intervalle de 


définition, est choisi a priori. Il n'y a donc que deux paramètres à ajuster : p et p. 
Les paramètres sont estimés par la méthode du maximum de vraisemblance : 


> (xj —x,)PLn(x; —X,) 
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La première équation est résolue par itérations successives, puis on résout la seconde. Pour cette 
distribution, les intervalles de confiances ne sont pas donnés. 


2.1.8 La loi de Poisson 


La mise en œuvre de l'ajustement à une loi de Poisson se distingue nettement de la procédure commune 
aux autres lois. 

Cette loi est utilisée principalement en hydrologie pour modéliser le nombre par an d'événements 
dépassant un seuil. Généralement les données se présentent donc sur la forme de deux colonnes : une colonne 
donnant les dates d'occurrence, la seconde les valeurs prises par la variable. Sur la figure 18, on trouve les 
données concernant les crues de la Garonne de 1913 à 1977, ayant un débit supérieur à 2500 m°/s. On se propose 
d'ajuster à une loi de Poisson le nombre de crue par an ayant dépassé ce seuil (4 en 1913, 3 en 1914, ... , 4 en 
1974, 1 en 1975, 2 en 1976 et 2 en 1977). Il suffit de sélectionner les années où se sont produites ces crues 
(fig.18) puis de lancer la commande "loi de Poisson" du menu "HYDROLAB"!. 
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Fig. 19 : introduction manuelle du nombre de modalité 


Fig. 18 : sélection des occurrences 


On obtient alors la feuille "POIS.XLS" dans laquelle vous devez introduire manuellement à la place du 
signe ?, le nombre de modalités possibles (voir figure 19). En effet, il n'est pas évident que le nombre de 
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modalités possibles soit 65 (1977-1913+1). Il se peut que la série présente des lacunes (nombre inférieur à 65), 
ou au contraire on peut imaginer que la série va au-delà de 1977 mais qu'il n'y a pas eu de dépassement du seuil 
(nombre supérieur à 65). Dés que l'on a rentré ce nombre de modalités, les feuilles de résultats se mettent à jour. 


Ajestement à une loi de Poisson 
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Fig. 20 : Feuille de résultats graphiques 
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Fig. 21 : Feuille de résultats numériques 


Sur le plan technique la loi de Poisson (ou loi des événements rares) indique que la probabilité de 
rencontrer k événements est donnée par : 


2.2 Régressions simples 


2.2.1 


Les variables à traiter sont obligatoirement situées dans deux colonnes, la première correspondant à la 


Sélection des données 


-p 
Proba(k) = e k. 


k 


k! 
u est le nombre moyen d'événement (ici 2,32=151/65) et on peut montrer que l'écart-type de k doit être 


égal à Ju (soit ici 1,52 = 42,32 , ce qui n'est effectivement pas éloigné de l'écart-type expérimental 1,63). 


variable à expliquer et la seconde à la variable explicative (si dans votre tableur cet ordre n'est pas 
respecté, il convient de le changer). Vous devez sélectionner uniquement les données numériques. La 
figure 22 vous donne trois exemples de sélection. 
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Fig. 22 : Exemples de sélection pour une régression simple 
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2.2.2 Exécution de la macro 


Bite de A ici TR émet Foie 1 


La figure 23 illustre le lancement d'une 
régression simple entre les pluies d'Antibes et 
Cannes sur la période commune 1959 - 1967. 

Une fois la macro lancé, vous trouverez 
j i e U na n u a deux feuilles analogues à celles décrites pour 
; : i les ajustements. 


g Une feuille intitulée "Regs.xlc" présente 
b: les points expérimentaux, la droite de 
Pist régression et l'intervalle de confiance sur 
1 l'estimation de la variable à expliquer 
= (première colonne et donc ici Antibes) en 
E7 fonction de la variable explicative (figure 24). 
j= Les valeurs numériques sont données 
r EN ES dans la feuille "Regs.xls" comme l'illustre la 
tal El) a $ figure 25. 
[bsi] Iim [E] Im Ig LIN peo ME TAN AL 


| Eemia iaaa | era | 
ner a rrn akon ce | LR 


Dar oi LS mure di | Ù à 


LME 


u 
A ed ne 
Fig. 24 : Graphique de la régression simple SEPP 


Fig. 25 : Valeurs numériques de la régression simple 


Soit Xk» Yk» kOx Et xOy les moyennes et écarts-types déterminés d'après les k observations communes ainsi 


que le coefficient de corrélation correspondant. 


Xps 


La régression de y en x s'écrit : 


Yx, représentant évidemment la moyenne conditionnelle des y liés par xj et si les hypothèses de normalité sont 
admissibles, y,, est la valeur la plus proche de y pour x; donné. 


Les bornes y, et y2 de l'intervalle de confiance à &% sur l'estimation de y à partir d'une réalisation x; de la 
variable explicative sont données par la relation suivante : 


p k-17 1 (x; -X Ÿ 
=y, tfh- EE hta ion 
yi2 = Fx, i P a k Loi(k—1) 


(formule tirée de YEVJEVICH 1982) 
Un tableau situé en bas de la feuille "Regs.xls" permet d'évaluer les valeurs de y connaissant x. 
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2.2.3 Personnalisation des résultats 


L= 


Kapara bre are ergh 


Cian de riha 
ag ime 


al | 
Fig. 26 : Personnalisation des résultats 

Comme le montre la figure 26, il suffit de reporter (taper au clavier ou coller) dans les cellules B11 et 
C11 les noms des variables à expliquer et explicative pour que toute l'interface utilisateur en tienne compte. Dans 
le tableau du bas, on a cherché à évaluer la pluie à Antibes pour l'année 1968 où il était tombé 1163 mm à 
Cannes et ce qu'il en serait à Antibes pour des pluies à Cannes de 500 et 1500 mm. On aurait pu également 
changer la borne de l'intervalle de confiance en remplaçant le 70% de la cellule F10 pour toute autre valeur qu'on 
aurait souhaitée. 


2.3 Régression double 


2.3.1 Sélection des données 


Les données sont disposées comme pour la régression simple mais sur trois colonnes ; la première 
colonne contenant toujours la variable à expliquer. 
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Fig.27 : Sélection des données et lancement d'une régression double 
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2.3.2 Exécution de la macro 
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Fig. 28 : Résultats numériques d'une régression double 


A cette feuille de calcul sont associées trois feuilles graphiques qui présentent les régressions partielles 
de la variable à expliquer par rapport aux deux variables explicatives ainsi que les reconstitutions de la variable à 
expliquer en regard de celles réellement observées (voir figure 30). 

Pour les aspects théoriques on se reportera à ce qui est dit au paragraphe 2.6 pour ce qui concerne les 
régressions multiples en général. 


2.3.3 Personnalisation des résultats 

Tout comme pour les régressions simples on peut personnaliser les feuilles en tapant les noms des 
variables (B8 à D8), en changeant les intervalles de confiance (122), en introduisant des valeurs des variables 
explicatives (F25 et G25), en ajoutant du texte sur les graphiques... 


0673 
0718 
pos 1900 | P.Chns | 


Nombre d'observations : w p 


945 “D -0,01 *P.Clanst+ -356552%8 


Fig. 30 : Personnalisation des résultats de la régression double 
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Fig. 31 : Graphiques personnalisés de la régression double 


2.4  Détections d'anomalies 


2.4.1 Sélection des données 


La figure 32 illustre le mode de 
sélection des données qui est analogue 
à celui de la régression simple. Par 
contre la convention est que se trouve 
dans la première colonne la variable à 
tester et que la seconde colonne 
contient la variable de référence 
réputée sans erreurs. 


Fig. 32 : Sélection des séries pour une détection des anomalies 
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2.4.2 Exécution de la macro 
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Fig. 33 : Résultats de la détection d'anomalies 


Cette macro a pour objectif de déterminer deux types d'anomalies : les anomalies ponctuelles et les 
anomalies systématiques. 
Dans le premier cas, on signale dans la colonne de gauche, les valeurs dont les résidus de régression 
dépassent une certaine probabilité d'être due au hasard. 
Si comme précédemment, les variables x et y sont corrélées et que les hypothèses d'une distribution de 
Gauss à deux dimensions sont acceptables, on aura les relations suivantes : 
y, = ax; +b et Yi—Yx =E 


i 


e; est une variable gaussienne de moyenne nulle et d'écart-type ©, V1- o? . Pour chaque valeur y;, on 


calculera le €; correspondant et sa fréquence théorique. Les valeurs de g; ayant des fréquences très rares 


correspondent à des y; douteux. On peut, par cette méthode, détecter des erreurs accidentelles qui n'apparaissent 
pas à l'étude des distributions marginales. Dans l'exemple de la figure 33 aucune valeur ne paraît anormale au 
seuil standard de 98%. 

Les anomalies systématiques sont détectées sur la base de l'analyse du cumul des résidus de régression : 
Comme on vient de le dire au paragraphe précédent, le résidu €; est une variable aléatoire gaussienne de 


moyenne nulle et d'écart-type GyV 1-6? . On définit alors la variable S;, cumul des i premiers résidus : 
i 
j=1 
E E 7 iln —1) 
Cette variable S; est une variable aléatoire de moyenne nulle et d'écart-type : 67 = O, Pg 
i i= 


Si on se fixe par exemple un intervalle de confiance à 98%, il y à deux chances sur cent pour que S; soit 
extérieur au segment : 


|- U09907, :+U0.9907, ] Soit encore : |- 2.3307, ,+2.3307, ] 


Lorsque j varie, le lieu des limites des segments ayant la même probabilité de contenir sont des ellipses 
passant par Z=0 pour i=0 et Z=0 pour i=n. 

Sur le graphique ci-dessous, on a porté les cumuls des résidus en fonction de i et la courbe restant à 
l'intérieur de l'ellipse à 98%, les deux séries paraissent homogènes. 
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2.4.3 Personnalisation des résultats 
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Sur la figure 34 on voit que après avoir 
collé les dates à partir de la cellule G12, 
indiqué les noms des variables (en C12 et 
D12) et modifié l'intervalle de confiance en le 
passant à 95% (cellule E9), la méthode 
détecte une anomalie en 1966 où la pluie à 
Antibes est de 1081.5 mm alors que compte 
tenu des 1166 mm à Cannes on aurait pu 
s'attendre à 949 mm. Cet écart de 132 mm a 
moins de 5% de chance d'être du au hasard et 
la macro nous indique que cette valeur est à 
vérifier. 


Sur la figure 35 nous avons traité toutes les 
pluies communes aux stations d'Antibes et 
Cannes (1959 - 1992 et à l'exception des deux -gi 
années 68 et 87) ; on remarque très nettement 
une rupture de stationnarité de par et d'autre des 
années 70 à 72. 


Dunectinn des anomaloa sromasiciy aee de Andes pas nipper à Canas 


En fait, il s'agit ici d'un déplacement de la 
station de (Cannes qui n'est donc pas la 
référence. 


heips agaat RES ce cha de contra I cum dhs ata 
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Fig. 35 : Détection d'une rupture de stationnarité. 


2.5 Analyse en composantes principales 


2.5.1 Préparation des données 
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Fig. 36 : Ouverture de la feuille modèle pour ACP 


La nature des données étant beaucoup 
plus complexe nous avons préféré passer 
par l'intermédiaire d'une feuille de calcul 
spécifique notée "Modacp.xls". Cette feuille 
est un document modèle dans lequel vous 
aller pouvoir regrouper tranquillement par 
le biais de "Copier - Coller" les 
informations à traiter. 


La première étape consiste à ouvrir 
dans le répertoire HYDROLAB le 
document "Modacp.xls" comme indiqué à 
la figure 36 


Cette feuille comporte un certain nombre 4 A B|CIDIEIFI6| M | 1 La eq 1M 
de cellules prévues pour l'introduction des E] oys F 
données. La ligne six, à partir de la |3 |  (dontnonmbes de vasisbtes sup) 0 pus 
3 ; FA (dont noeitee d'observations sup) Ò «uplus 
colonne B contiendra le nom des variables O emptus 
(des alphanumériques ou des nombres) de 
même la colonne A à partir de la ligne sept 
contiendra le nom des observations (des 
alphanumériques ou des nombres). La 
plage des cellules à gauche et en dessous 
de la cellule B7 est destinée à recevoir les 
données numériques (il ne faut absolument 
pas qu'il y ait de lacunes: blancs ou Fig. : 37 Aspect de la feuille "Modacp.xls" à l'ouverture 


alphanumériques). 
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Fig. : 38 Sélection des données pour l'ACP j 
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La figure 38 donne un exemple de sélection des données à traiter par ACP, ici, on a sélectionné 8 variables (les 
postes pluviométriques et 13 observations (les pluies annuelles de 69 à 81). Une fois cette plage de cellules 
sélectionnée, on vient de la copier puis de la coller dans la feuille "Modacp.xls" comme l'illustre la figure 39. 


Microsoft Excel - Modacp xls 
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Fig. :39 Collage des données d'ACP 


Il est de loin préferable de procéder par un collage spécial des valeurs, en effet certaines de vos données 
de départ peuvent être des résultats de formules. Dans l'exemple ci-dessus le transfert a été fait en une seule fois 
(nom des variables, nom des observations et les valeurs numériques), mais vous pouvez procéder en plusieurs 
étapes si cela vous est plus commode. Ces différents collages effectués vous devriez obtenir quelque chose 


d'analogue à la figure 40. 
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Fig. : 40 Introduction des paramètres de l'ACP 


Le nombre de variables supplémentaires (ce sont 
des variables qui n'entreront pas dans le calcul des 
composantes mais qui seront néanmoins projetés 
dessus). Ces variables supplémentaires 
correspondent aux colonnes les plus à droite. 

Le nombre d'observations supplémentaires (ce sont 
des observations qui n'entreront pas dans le calcul 
de la matrice de coefficients de corrélation mais 
qui seront néanmoins projetés sur les 
composantes). Ce seront les observations figurant 
dans les dernières lignes. 

Enfin, indiquez le nombre de composantes désirées 
; ce nombre doit être notablement inférieur au 
nombre de variables entrant dans les calculs 
(nombre total des variables moins nombre de 
variables supplémentaires). 


2.5.2 Exécution de l'ACP 


LE Ne vous inquiétez pas du format des 
colonnes (des nombres remplacés par des #, 
des noms de variables illisibles), qui sera 
géré automatiquement. Par contre il vous 
reste à saisir : 

e le nombre total de variables à traiter, 
(dans la feuille on vous indique la 
valeur maximale compte tenu des 
données introduites, mais vous pourriez 
ne vouloir traiter que les x<ë8 
premières). 

e Le nombre total d'observations (dans la 
feuille on vous indique la valeur 
maximale compte tenu des données 
introduites, mais vous pourriez ne 
vouloir traiter que les y<13 
premières). 
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Fig. :41 Lancement de l'ACP 


Il suffit alors de lancer la commande ACP comme l'indique la figure 41. 
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po | 


Fig. :42 Boites de dialogue de l'ACP 


À la première boite de dialogue il convient de 
répondre oui. (la macro utilise un fichier temporaire 
déjà existant "Donacp.txt". 


Pour la seconde boite de dialogue il convient de 
répondre non, en effet, HYDROLAB est conçu pour 
tourner indifféremment dans les versions 4, 5 et 97 
d'Excel, les fichiers sont donc tous à sauvegarder en 
version 4 (compatibilité ascendante). 
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Vous verrez apparaître alors la fenêtre DOS du 
module exécutable. Cette fenêtre permet de suivre 
sommairement le déroulement des opérations. Il se peut 
que l'on reste bloqué dans cette fenêtre si l'on a commis 
des erreurs dans la construction du tableau de données 
(Modacp.xls). Les erreurs les plus banales sont un 
tableau comportant des valeurs non numériques et des 
variables étant des combinaisons linéaires des autres. 
En cas de blocage, revenir à la saisie des données. 

Suivant les versions, cette fenêtre se ferme 
automatiquement ou sinon il vous suffit de cliquer sur 
la croix en haut à droite de la fenêtre (figure 43) 


2.5.3 Visualisation des résultats de l'ACP 
B erher piton fichage lréerten Format Quik |Hygroib Qornées Faute t 
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2 Code commune 4 fa] n Loic Poison 
4 Code beute l | l LoiRscre-rormels 
x pn INA x 
6 i t LH ws Lawetut 
BE 4 1 
gM M Psmemens _ 
ts x nm us 
10 tx i 15 
11 ve pe 10 
12 1%2 ét M6 La Reg eerion double 
Fig. : 44 Visualisation des résultats de l'ACP 
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Fig. : 43 Fenêtre du module exécutable de l'ACP 


Il vous suffit alors de consulter les résultats. 
Pendant quelques instant la commande va récupérer 
les résultats du module exécutable et les remettre sous 
la forme d'un tableau Excel nommé "Resacp.xls". Ce 
fichier sera écrasé à chaque nouvelle utilisation. Une 
boite de dialogue vous demandera l'autorisation de 
remplacer l'ancienne version, il convient de répondre 
par oui. 
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Fig. : 45 Exemple de résultats d'une analyse en composantes principales 

La programmation de l'analyse en composantes principales s'inspire des ouvrages de LEBART & al. 
(1979) et de LEGRAS (1972). Ce programme a été écrit en Basic compilé, seul le module exécutable est 
indispensable mais vous pouvez consulter la source "Acpjpl.bas" dans le répertoire HYDROLAB. 

Il faut noter que nous avons choisi de présenter les projections des observations sur les composantes de 
façon centrée réduite. Dans certains autres programmes ces projections sont données avec une variance égale à la 
valeur propre correspondante. Il ne faut donc pas s'alarmer si vous constatez des écarts avec les résultats d'autres 
programmes ; la convention n'est pas la même et il suffit d'un terme multiplicatif pour homogénéiser les 


résultats. 


Si vous voulez conserver ces résultats il suffit de les enregistrer sous un autre nom. 
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2.6  Régressions multiples 


2.6.1 Préparation des données 
Ouvre HE 
Page dansi [2 Ho ou =] ©) aaa = m e La préparation des données est tout à 
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fait analogue à ce que nous avons décrit 
pour l'analyse en composantes principales. 
r Nous passons par l'intermédiaire d'une 
feuille de calcul spécifique notée 
"Modreg.xls". Cette feuille est un document 
modèle dans lequel vous aller pouvoir 
regrouper tranquillement par le biais de 
zj "Copier - Coller" les informations à traiter. 


La première étape consiste à ouvrir 
dans le répertoire HYDROLAB le 
document "Modreg.xls" comme indiqué à la 


Fig. 46 : Ouverture de la feuille modèle pour régression multiple figure 46 


Cette feuille comporte un certain nombre de 
cellules prévues pour l'introduction des données. 
La ligne trois, à partir de la colonne B 
contiendra le nom des variables (des 
alphanumériques ou des nombres) de même la 
colonne A à partir de la ligne quatre contiendra 
le nom des observations (des alphanumériques 
ou des nombres). La plage des cellules à gauche 
et en dessous de la cellule B4 est destinée à 
recevoir les données numériques (il ne faut 
absolument pas qu'il y ait de lacunes: blancs ou 
alphanumériques). Il est indispensable de mettre 
dans la première colonne la variable à expliquer, 
puis les variables explicatives dans les colonnes 
suivantes. 


1 y 

2 Lieu-dit Place Neuve 

3 Code commune Sà 4 
4 Code lieu-dit 1 2 
5 jx 18770 18522 
6 ly 9768 9878 
Ag 1000 74 
8 Année 33 32 
9 1959 2034 963 
10 1960 2324 1127 
11 1961 1165 b44 
12 1962 1193 621 
13 1963 1668 638 
14 1964 132 658 
15 1965 1062 592 
16 1966 1276 1082 
17 1967 74 358 


mA? qar 671 


A] Modreg.xis 
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Fig. : 47 Aspect de la feuille "Modreg.xls" à l'ouverture 
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Fig. : 48 Sélection des données pour la régression multiple 
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La figure 48 donne un exemple de sélection des données où on cherche à expliquer les pluies à Coursegoule (la 
colonne correspondante de la feuille de données "pluies06.xls" a été permutée) à partir de sept autres stations du 
département. La régression est étudiée à partir de 13 observations (les pluies annuelles de 69 à 81). Une fois cette 
plage de cellules sélectionnée, on vient de la copier puis de la coller dans la feuille "Modreg.xls" comme l'illustre 


la figure 49. 
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Fig. :49 Collage des données de la régression multiple 


Il est de loin préférable de procéder par un collage spécial des valeurs, en effet certaines de vos données 
de départ peuvent être des résultats de formules. Dans l'exemple ci-dessus le transfert a été fait en une seule fois 
(nom des variables, nom des observations et les valeurs numériques), mais vous pouvez procéder en plusieurs 
étapes si cela vous est plus commode. Ces différents collages effectués vous devriez obtenir quelque chose 


d'analogue à la figure 50. 


Fig. : 50 Introduction des paramètres de la régression multiple 
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Ne vous inquiétez pas du format des 
colonnes (des nombres remplacés par des #, 
des noms de variables illisibles), qui sera 
géré automatiquement. Par contre il vous 
reste à saisir : 


le nombre total de variables, (le 
nombre de variables explicatives plus 
un). Si vous tapez un nombre n 
inférieur à celui proposé (par exemple 
six au lieu de huit), la première 
variable sera expliquée par les n-1 
variables suivantes. 
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:51 Lancement de la régression multiple 
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2.6.2 Exécution de la régression multiple 


Il suffit alors de lancer la commande Régression Multiple comme l'indique la figure 51. 


à) Un Se rom C HIDAQLATONRES TIT ete dé à cat emplacement Voleros le 
d mga? 


Caed] 
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Fig. :52 Boites de dialogue de la régression multiple 


Vous verrez apparaître alors la fenêtre DOS du 
module exécutable. Cette fenêtre permet de suivre 
sommairement le déroulement des opérations. Il se peut 
que l'on reste bloqué dans cette fenêtre si l'on a commis 
des erreurs dans la construction du tableau de données 
(Modreg.xls). Les erreurs les plus banales sont un 
tableau comportant des valeurs non numériques et des 
variables étant des combinaisons linéaires des autres. 
En cas de blocage, revenir à la saisie des données. 

Suivant les versions, cette fenêtre se ferme 
automatiquement ou sinon il vous suffit de cliquer sur 
la croix en haut à droite de la fenêtre (figure 53) 
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Fig. : 54 Visualisation des résultats 
de la régression multiple 
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A la première boite de dialogue il convient de 
répondre oui. (la macro utilise un fichier temporaire 
déjà existant "Donreg.txt". 


Pour la seconde boite de dialogue il convient de 
répondre non, en effet, HYDROLAB est conçu pour 
tourner indifféremment dans les versions 4, 5 et 97 
d'Excel, les fichiers sont donc tous à sauvegarder en 
version 4 (compatibilité ascendante). 
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Fig. : 53 Fenêtre du module exécutable de la 
régression multiple 


Visualisation des résultats de la régression multiple 


Il vous suffit alors de consulter les résultats. 
Pendant quelques instant la commande va récupérer 
les résultats du module exécutable et les remettre sous 
la forme d'un tableau Excel nommé "Resreg.xls". Ce 
fichier sera écrasé à chaque nouvelle utilisation. Une 
boite de dialogue vous demandera l'autorisation de 
remplacer l'ancienne version, il convient de répondre 
par oui. 
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Fig. : 55 Exemple de résultats d'une régression multiple 
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La programmation de la régression multiple s'inspire des travaux de TOUCHEBOEUF de LUSIGNY . Ce 
programme a été écrit en Basic compilé, seul le module exécutable est indispensable mais vous pouvez consulter 
la source "Regjpl.bas" dans le répertoire HYDROLAB. 

Nous avons adopté la démarche dite régressive (backward elimination), elle consiste dans un premier 
temps à faire une régression multiple sur l'ensemble des variables explicatives envisagées. Une fois ce premier 
résultat établi, on vérifie la signification des coefficients de corrélation partielle (en fait, il suffit de tester le plus 
petit coefficient de corrélation partielle en valeur absolue). On teste cette signification par le test de Student au 
seuil de 5% (seuil choisi à priori est non modifiable). Ensuite on teste le coefficient de corrélation multiple par le 
test de Fisher-Snedecor au seuil de 5%. Si l'un de ces deux tests est négatif on supprime la variable explicative 
correspondant au plus petit coefficient de corrélation partielle, et on recommence avec une variable explicative 
en moins. Le programme s'arrête lorsque tous les tests sont positifs. 

Sur la figure 55, vous trouvez tous les résultats intermédiaires justifiant la démarche. Successivement les 
variables explicatives Cannes, Colomar, la Colle, Contes et enfin Entraunes sont éliminées. Dans ce cas de 
figure, la régression finale revient à expliquer les pluies à Coursegoule par les seules observations d'Antibes et de 
Clans. La régression à pour expression : 

Pcoursegoules = 0.86 PA ntibes + 0.67Pcrans + 65.85 


Le coefficient de corrélation multiple est de 0.96 et les coefficients de corrélation partielle avec Antibes et 
Clans sont respectivement 0.80 et 0.78. 
Si vous voulez conserver ces résultats il suffit de les enregistrer sous un autre nom. 


2.7 Comblement de lacunes (Bouche trous) 


L'objectif de cette commande est de combler des lacunes dans une matrice de données. Les variables sont 
comme toujours disposées suivant des colonnes et les observations suivant des lignes. Dans cette matrice de 
données certaines valeurs absentes seront remplacées par un texte contenant la lettre a (par exemple : "Abs.", 
"Lacune", ...). On trouvera plus loin la justification du calcul il suffit pour l'instant de savoir que le nombre de 
composantes principales à prendre en compte est au libre choix de l'utilisateur. Pour des données 
pluviométriques par exemple, il est recommandé de ne retenir que les composantes présentants une structure 
spatiale. L'opérateur a également le choix du nombre d'itérations à effectuer (10 à 15 sont généralement 
suffisantes). Pour plus de détail sur cette méthode on peut se référer à l'extrait du cours LABORDE (1998) donné 


en annexe I. 


2.7.1 Préparation des données 
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= H 3) 3 n D on o e S n les noms des observations. La figure 56 


illustre une telle sélection de huit variables et 
douze observations. Cette plage étant 
sélectionnée, vous allez la copier et la coller 
ié dans la feuille "Modbt.xls" qui se trouve dans 
714 le répertoire HYDROLAB. Une fois ce 
collage fait, il vous suffit de compléter les 
1155 quatre premières lignes du tableau comme le 
montrent les figures 57 et 58. 
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Fig. 57 : Ouverture de la feuille modèle "Modbt.xls" Fig. 58 : La feuille modèle prête à être utilisée 
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2.7.2 Exécution de Bouche trous 


Lorsque la feuille modèle est complétée telle qu'à la figure 58, il suffit de lancer la commande "Bouche trous" du 


menu déroulant HYDROLAB. 
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Fig. 59 : Lancement de "bouche trous" 


Vous verrez apparaître alors la fenêtre 
DOS du module exécutable. Cette fenêtre 
permet de suivre sommairement le déroulement 
des opérations. Il se peut que l'on reste bloqué 
dans cette fenêtre si l'on a commis des erreurs 
dans la construction du tableau de données 
(Modbt.xls). Les erreurs les plus banales sont un 
tableau comportant des valeurs non numériques 
sur toute une ligne ou toute une colonne et des 
variables étant des combinaisons linéaires des 
autres. En cas de blocage, revenir à la saisie des 
données. 


Fig. 60 : Fenêtre du module exécutable de "Bouche trous" 


Une fois lancée la commande, deux boites de dialogue 


Vous seront soumises : 


A la première boite de dialogue il convient de 


répondre oui. (la macro utilise un fichier temporaire déjà 
existant "Donbt.txt". 


Pour la seconde boite de dialogue il convient de 


répondre non, en effet, HYDROLAB est conçu pour 
tourner indifféremment dans les versions 4, 5 et 97 
d'Excel, les fichiers sont donc tous à sauvegarder en 
version 4 (compatibilité ascendante) 


Suivant les versions, cette fenêtre se ferme automatiquement ou sinon il vous suffit de cliquer sur la 


croix en haut à droite de la fenêtre (figure 60) 


2.7.3 Visualisation des résultats de Bouche trous 


cher gen ancrage puetien Format Qui [Haroub Qornées faite 2 
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Une fois la commande "Bouche trous” exécutée, 
vous pourrez visualiser les résultats en lançant la 
commande "Lecture des résultats de Bouche trous" du 
menu déroulant d'HYDROLAB. 

Pendant quelques instant la commande va 
récupérer les résultats du module exécutable et les 
remettre sous la forme de deux tableaux Excel nommés 
"DetResbt.xls" et "Resbt.xls". Ces fichiers seront 
écrasés à chaque nouvelle utilisation. Une boite de 
dialogue vous demandera l'autorisation de remplacer 
l'ancienne version, il convient de répondre par oui. 

Les résultats sont représentés d'une part dans un 
fichier contenant la matrice originale (Resbt.xls) mais 
où les lacunes sont remplacées par leurs estimations et 
d'autre part dans un fichier donnant le détail des 
évolutions des estimations d'une itération à l'autre 
(Detresbt.xls). Ce fichier ne sert qu'à vérifier la 
convergence du processus. 
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2.8 Variographie 


La mise au point de cette commande s'explique par l'utilisation fréquente du logiciel SURFER dans les 
synthèses cartographiques. Ce logiciel est interfacé avec Excel, mais il ne comprend pas de module d'analyse 
géostatistique. L'objectif de cette commande est donc d'effectuer une modélisation des variogrammes 


expérimentaux de façon à guider le choix dans les paramètres de krigeage de SURFER. Ce texte s'inspire 
principalement des travaux de C. OBLED 1986. 


2.8.1 Sélection des données 


Chaque point de mesure est caractérisé par ses coordonnées cartésiennes et par la valeur mesurée en ce 
point de la fonction dont on veut étudier la fonction structure. Il convient donc de sélectionner dans une feuille 
Excel les coordonnées ainsi que les mesures sous la forme de trois colonnes dont la première contient les 
abscisses, la deuxième les ordonnées et la troisième la variable. Dans la première ligne on sélectionnera les noms 
de ces colonnes. La figure 62 illustre la sélection des coordonnées des postes pluviométriques des Alpes 


maritimes disposant d'au moins trente années de mesures, et de la précipitation moyenne interannuelle. 
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Fig. 63 : Lancement de la commande " Calcul des 
construction du variogramme accroissements pour le variogramme 


Fig. 62 : sélection des données pour la 


2.8.2 Constitution des couples 


Cette sélection effectuée, on lance la commande "calcul des accroissements pour le variogramme" 
comme le montre la figure 63. Cette commande permet de constituer la feuille "Vario.xls" dans laquelle sont 
récapitulés tous les couples de points de mesure envisageables avec les distances entre points, les azimuts relatifs 
et les accroissements des mesures. 

Comme le montre la figure 64, en 
colonne E et F sont indiqués les principales 
caractéristiques de l'échantillon à savoir sa 
taille, la moyenne, la variance et l'écart-type 
des valeurs, le nombre de couples possible 
et la distance maximale entre deux points de 
mesure. Dans les colonnes AA et AB, vous 
avez la possibilité de jouer sur les 
paramètres déterminants les couples à 
prendre en compte dans la constitution du 
variogramme expérimental. 

La première précaution est de prendre 
maximale d'analyse largement inférieure à 
la distance maximale (la moitié, voir le 


Van vd de 


| Dasare mas arabe] 1751} 


Lout-rael 112 
kheo 


NS RONSE EST a Tarn 
$ 
3 


=i tiers). Ayant choisi cette distance, vous 
a 23| 14718) 10014 e2 pouvez également fixer le pas de distances 
ET an| verre] vai) en séparant les couples en classes. Enfin, il est 
D va me, vo possible de ne prendre en compte que les 


couples tels que leur azimut relatif soit 


Fig. 64 : Description de la feuille "Vario.xls" compris entes 
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Fig. 65 : Paramètres d'un variogramme anisotrope 


2.8.3 Calcul du variogramme expérimental 


A7, - tolérance et AZ, + tolérance 

Ainsi que l'indique la figure 65, il est 
donc possible de construire des 
variogrammes expérimentaux anisotropes. 
Attention pour l'azimut, la convention est de 
donner l'angle en degré, à partir de l'axe des 
abscisses et en tournant dans le sens 
trigonométrique (sur la figure on a 
sensiblement Azo=30° et une tolérance de 
15°). 

À priori, les paramètres proposés sont 
ceux d'un variogramme omnidirectionnel 
(tolérance de 90°) calculé sur la moitié du 
champ. Le pas d'exploration des distances 
est le 20°* de la moitié du champ (20 
classes de distances). 


Fig. 66 : Choix des caractéristiques du variogramme. 


Une fois choisie les paramètres du variogramme (ici 
Fig. 66, on a pris une tolérance de plus ou moins 45° 
autour d'une direction moyenne à 10°, la distance 


maximale d'analyse est de 400 et l'on a des classes de 
40), on peut lancer la construction des variogrammes 
théoriques 


expérimentaux et comme 


la fig. 67. 


l'indique 


[rm rare sc #0 r 
Air 


Fig. 68 Feuille récapitulative des caractéristiques des variogrammes 


. 67 :Lancement du calcul des variogrammes 


théoriques et expérimental. 


(1 Nanas a LORS ETES 


Fig. 69 Représentation antique d des 
variogrammes théoriques et expérimental 


Le lancement de cette commande permet la construction de deux feuilles. La première comme l'illustre 
la figure 68 récapitule les différents éléments constitutifs du variogramme expérimental (pour chaque catégorie 
on trouve successivement l'effectif, la distance moyenne et le y(h) expérimental), et dans la partie droite, les 
caractéristiques des variogrammes théoriques ( type, portée palier et pépite), la colonne "erreur type" caractérise 
les écarts entre le variogramme expérimental et les théoriques. Cette colonne est un guide pour l'optimisation des 
paramètres du variogramme théorique à ajuster. La seconde feuille Fig. 69, présente les résultats précédents sous 


forme graphique. 
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2.8.4 Ajustement du variogramme théorique 


MOT DU MR UT Had PR RE Avant d'ajuster les paramètres il convient 
tout d'abord de choisir un des quatre modèles 
théoriques proposés. Dans l'exemple de la fig. 
69, il paraît raisonnable de retenir un modèle 
gaussien (l'effectif de la première catégorie 
n'est que de 1 mais la suite semble indiquer un 
comportement très lisse à l'origine). Il nous 
reste alors que trois paramètres à ajuster. La 
première méthode consiste à modifier 
manuellement les paramètres et à tâtonner au 
regard des graphiques. 

La seconde méthode consiste à utiliser le 
Solveur d'Excel comme l'illustre la Fig. 70. 

On obtient alors une boite de dialogue dans 
laquelle on va préciser ce que l'on veut 
optimiser, les paramètres susceptibles d'être 
modifiés et d'éventuelles contraintes. 


r Non amp 


La fig. 71 montre l'ajustement à un 
modèle gaussien obtenu par minimisation 
de l'erreur type correspondante (cellule 
US) en faisant varier les portée, palier et 
pépite (cellules RS à T5) sous la seule 
contrainte que la pépite soit positive 
(cellule T5). CAES TO TT -an 
L'optimisation est effectuée en cliquant sur [a He me œur 
le bouton "résoudre". ‘| == 


TEETE 0 a 


Fig. 71 Choix des paramètres du Solveur. 
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Fig. 72 Résultat de l'ajustement par le Solveur. 


Dans la mesure où l'on avait choisi le modèle gaussien, nous avons allégé le graphique de la Fig. 72 en 
supprimant les autres modèles. Il suffit pour cela de cliquer sur la courbe correspondante et de taper la touche 
"supprimer". 


2.8.5 Construction d'un nouveau variogramme 


Si sur les même données, on désire construire un nouveau variogramme expérimental (changement de 
direction, de tolérance, de pas de distance...) il convient de fermer les deux feuilles "variomod.xls" et 
“variomod.xlc". On revient alors à la feuille "vario.xls" dans laquelle on indiquera les nouveaux paramètres de 
construction du variogramme expérimental puis on procédera à nouveau comme indiqué au paragraphe 2.7.3. 


2.8.6 Définitions et rappels théoriques 
Pour plus de renseignements on peut se référer à l'article de C. OBLED (1986). 
Les quatre modèles proposés correspondent aux définitions suivantes : 


E Pali 
Modèle linéaire : y(h) = Pépite +h = = 
o 


2 


tée 
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3 
y(h) = Pépite + Palier | — LS 


Si h{ Portée 
2 Portée 2 E ( 


Modèle Sphérique : 
y(h)= Pépite + Palier Si h) Portée 


3h 
Modèle Exponentiel : y(h) = Pépite + pated —e Portée | 


{ 1.732 h | 
Modèle gaussien : ylh) = Pépite + Palier| 1 — e Portée 


Comme le montre la Fig. 73, les définitions 
des variogrammes théoriques sont telles que 
lorsque la distance h est égale à la portée, le 
terme y(h) est sensiblement égal à 
Pépite + Palier. 

Ce que nous appelons "portée" est donc ce 
que d'autres appelle parfois "portée pratique" ou 
encore "distance de décorrélation". 

Ce choix paraîtra peut être curieux à des 
puristes mais il est justifié par son côté pratique. 
En effet, la "portée" que vous ajusterez sous 
HYDROLAB, correspond au paramètre 
"Range" nécessaire au krigeage sous SURFER. 


Fig. 73 : Variogrammes théoriques utilisés 


3 Les fonctions personnalisées 


Lo RaDio Soar 
AYDROLAB, xiaiCoef. _de_00 : . r 
YDROLAR ia Durée pu — Ces fonctions sont accessibles par le menu "Insertion" 
VIROLAR aF. do Gass 


NORCLAS iiair Aa N Ferrer Gnidia où en choisissant "Fonction" on accède à la boite de 
ru eh dialogue de la Fig. 74. En choisissant dans les catégories 
Rd de fonctions celles nommées "Personnalisées", on a accès 
ape zi en colonne de gauche à différentes fonctions 
nn seu re d'HYDROLAB. Elles sont rangées par ordre alphabétique 


D o e n IEE et précédées du préfixe "HYDROLAB.xla!". 


e) Ca y mu | 


Fig. 74 : Appel des fonctions personnalisées 
3.1 Fonctions statistiques 


3.1.1 U._de_Gauss 


Cette fonction permet d'évaluer la valeur de la variable réduite de Gauss U correspondant à une fréquence au 
non dépassement F. Elle ne possède qu'un seul argument d'entrée F qui doit être compris entre 0 et 1. 
Pour les justifications on peut se reporter au paragraphe 2.1.2. 
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3.1.2 


=HYDROLAB.xalU._de_Gauss(A2) 
0,02 =HYDROLAB xislU._de_Gauss(A) 
0,05 l=HYDROLAB xalU._de_Gauss(A4) 
01 (=HYDROLAB xalU._de_Gauss(A5) 
02 |=HYDROLAB xalU._de_Gauss(Af) 
05 =HYDROLAB xalU,_de_Gauss(AT) 
08 |=HYDROLAB xalU._de_Gauss(A8) 
09  =HYDROLAB.1alU,_ de Gauss(A9) 
095 =HYDROLABxdalU._de_Gauss(A10) 
038 =HYDROLAB xialU._de_Gauss(A!l) 
099 =HYDROLAB xialU._de_Gauss(A12) 


Fig. 75 : Exemple d'utilisation de la fonction U_de_Gauss 
(à gauche affichage des fonctions, à droite affichage des valeurs numériques) 


F._de_Gauss 


-2,326785333 
-2,054188589 


-1,64521144 


-1,281728757 
-0,841456717 


1,01007E-07 
0,841456717 
1,281728757 

1,64521144 
2054188589 
2,326785333 


Cette fonction est l'inverse de la précédente où connaissant une valeur de la variable réduite U, on peut obtenir 
une estimation de sa fréquence au non dépassement F. 
Pour les justifications on peut se reporter au paragraphe 2.1.2. 


U 


-2,32678533255897  (=HYDROLAE xlalF._de_Gauss(A2) 
-2,0541885887219 =HYDROLAB xlalF._de_Gauss(A3) 
-1,64521144014382  =HYDROLAB-xialF._de_Gauss(A4) 
-1,28172875650271  =HYDROLAB xlalF._de_Gauss(A5) 
-0,841456717354784  =HYDROLAB.slalF._de_Gauss(A6) 
L01006675468085E-07 =HYDROLAB.xlalF._de_Gauss(A7) 
0,841456717354784  =HYDROLAB xlalF._ de Gauss(A8) 
128172875650271  =HYDROLAB xaiF._de_Gauss(AS) 
L64521144014381  =HYDROLAB xlalF._de_Gauss(Al0) 
2,0541885887219 =HYDROLAE xlalF._de_Gauss(Al 1) 
232678533255897  =HYDROLAB xalF._de_Gauss(A12) 


3.1.3 


Cette fonction donne la fréquence au non dépassement d'une valeur F de Fisher Snedecor pour des degrés de 


-2,326785333| 
-2,054188589| 
-1,64521144| 
-1,281728757| 
-0,841456717| 
1,01007E-07 
0841456717. 
1281728757) 
164521 144) 
2054188589) 
2326785333) 


Fig. 76 : Exemple d'utilisation de la fonction F_de_Gauss 
(à gauche affichage des fonctions, à droite affichage des valeurs numériques) 


F. de F. de Fisher Snedecor 


liberté v1 et v2. 


La variable x = 


ee) 2 
9v, 9v; 

Eo ER 

9v; 9v, 


0,010011689 
0,019819783 
0,049778124 
0,100097926 
0.20014752 
0,50000004 
0,79985248 
0899302074 
0950221876 
0,980180217 
0989988311 


F de Fisher Snedecor. Pour plus de détail on se référera à ABRAMOWITZ (1964), page 947. 
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Lil =HTORDOLAB due de FFisotrér ideaa En BE BAG) no BA 0a 0m AN) AFA 
5 SHIDECLAE dut de F Fikh Arries ET BE EAT) D PS OSSI 0m aM WA 
X =HFOROLA B daf de FFi der U OF OI OA) Ml) MF 


AE 
Fig. 77 : Exemple d'utilisation de la fonction F. de F. 


de Fisher Snedecor 
(à gauche affichage des fonctions, à droite affichage des valeurs numériques) 


suit une loi de Gauss et c'est ainsi que l'on attribue une fréquence au 
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3.2 Fonctions liées à l'ETP et l'ETR 


3.2.1 Durée astronomique du jour (Duréejour) 


Cette fonction donne la durée du jour exprimée en heures en fonction de la date et de la latitude exprimée en 
degré et dixièmes (+ pour le nord et - pour le sud). 
| DORE) 


| A | B c | 
“Date Lantade H (heures) 


ÆUDa_e Latitude H (heures) 
21/0798 4271 1468048642 


2 135997 4271 =HYDROLAE da!Duréejour(B2,A2) 
2 21/07/98 4271! 9446538902 
18/01/99, 4271, 9437347098 


[8135997 4271 "HYDROLAB xdalDurésjour(B3,A3) 
18/01/99, -4271 14,81427023| 


36178 4271  =HYDROLAB xlalDuréejour(B4,A4) 

3617 4271 =HYDROLAE xdalDuréeyour(B5,A5) 
Fig. 78 : Exemple d'utilisation de la fonction Duréejour 

(à gauche affichage des fonctions, à droite affichage des valeurs numériques) 


Pour plus de détail on se reportera à l'annexe IT présentant un extrait de LABORDE (1997). 


3.2.2 Radiation globale théorique (Iga) 


Cette fonction donne la radiation globale théorique (en calories par cm? et par jour) en fonction de la date et de 
la latitude exprimée en degré et dixièmes (+ pour le nord et - pour le sud). 


A | 6 | c S NE SE E 
M [Date Latitude Iga M iDate Latitude Iga 
235796 45 =HYDROLAB. xla!Iga(B2.A2) 45, 261,5404155 
235855 45 =HYDROLAB xlallga(B3,43) 45. 523496272 
4135916 45 =HYDROLAB xla!Iga(B4. A4) 45, 898,0937273 
SM 35977 45 =HYDROLAB xlallga(B5, AS) 45 10135631 
6 (36039 45  =HYDROLAE xlallga(B6.A6) s Taua 
F (36100 45 =HYDROLAB xlallga(B7,A7) > 


Fig. 79 : Exemple d'utilisation de la fonction Duréejour 
(à gauche affichage des fonctions, à droite affichage des valeurs numériques) 


Pour plus de détail on se reportera à l'annexe I présentant un extrait de LABORDE (1997). 


3.2.3 ETP Penman originale (Penman) 


Cette fonction permet d'évaluer l'ETP Penman exprimée en mm/jour en fonction de la date, de la latitude (en 
degré et dixièmes), de l'humidité relative (en %), de l'insolation journalière (en heures et dixièmes), du vent 


moyen à 2 mètres (en m/s), de l'Albédo et de la température (en degré Celsius). 
ES EA A 28 e R ET A EO, EE en MAS r E A a a a: S 
1 
Dae Latitude Albédo temp. vent bumedté [nso] ETP Penman 


Date Lattade Albéde temp vem humidité insol ETP Penman 
liss 45 02 45 24 N 77 =HYDROLAB xalPereman(B3,A3,C3,D3%E3,03.F3) 2.2 


louise 45 GET 2 77 S 
W5? 45 02 23 52 59 129 =HYDROLAB. xda! Pernan(PA. A4.C4.D4.F4.C4,F4) 010698 4 52 S 129 7482 
S haos 45 02 BI 24 42 87 =HYDROLAB xalPerean(B5,A5.C5,D5,E5,05,FS) 010896 45 24 4a 8&7 6535 


Fig. 80 : Exemple d'utilisation de la fonction Penman 
(à gauche affichage des fonctions, à droite affichage des valeurs numériques) 


Pour plus de détail on se reportera à l'annexe I présentant un extrait de LABORDE (1997). 


3.2.4 ETP Penman modifiée FAO (FAO) 

Cette fonction permet d'évaluer l'ETP Penman modifiée par la FAO exprimée en mm/jour en fonction de la date, 
de la latitude (en degré et dixièmes), de l'humidité relative (en %), de l'insolation journalière (en heures et 
dixièmes), du vent moyen à 2 mètres (en m/s), de l'Albédo et de la température (en degré Celsius). 


EST B E AN E G H ] Ao] 
A Altitude = 755 Dlatmoie= 755 

P Latitude Albédo temp. vent Inerudté inso FAO E Date Latitute Albédo temp. vent humiisė iasol FAO 
O EL 2 45 24 M 77 =HYDROLAB mal FACO(ERAJ;C}DI;E3;03;F3) Ba ons 45 42, 4,5 24 T 77 am 
A 35w47 45 0,2 213 52 59 129 =HYDROLAB adal FACNEA,A4,C4,D4,E4,C4,F4) 4! 0106/93 45 o 3, 52 59, 129! 8,144 
| 4 45 2 239 24 42 87 =HYDROLAE MalFACNBS ASCS.DS,ES.05,F5) Sn onas 45 24 a T 7,2% 


Fig. 81 : Exemple d'utilisation de la fonction FAO 
(à gauche affichage des fonctions, à droite affichage des valeurs numériques) 
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Pour plus de détail on se reportera à l'annexe II présentant un extrait de LABORDE (1997). 


3.2.5 ETP Penman modifiée Mc Culloch (McCulloch) 


Cette fonction permet d'évaluer l'ETP Penman modifiée par Mc Culloch exprimée en mm/jour en fonction de la 
date, de la latitude (en degré et dixièmes), de l'humidité relative (en %), de l'insolation journalière (en heures et 
dixièmes), du vent moyen à 2 mètres (en m/s), de l'Albédo, de l'altitude (en m) et de la température (en degré 
Celsius). 


75 
Due [CuB Abiio tenp vert Hassètié Itael ETP MeOdech 
DNUS 4S O2 45 24 72 17 OION 
DOSE 4S Q2 HI SEO 59 119) 443642 
owe 45 02 289 24 42 A7 7e 


Fig. 82 : Exemple d'utilisation de la fonction McCulloch 
(à gauche affichage des fonctions, à droite affichage des valeurs numériques) 


357% 45 02 4,5 A R 17 ÆHYDROLAB shiMeOeMBSASCI,D3,E301,F38E81) 
35347 45 a2 21,3 152 99 129 ÆHYDROLAS slalMeCuloeh{B4,A4,C4,D4,E4,04,F4.$B81) 
ue 45 02 2395 124 42 47 =HYDROLAB shalMeOeb{BS,AS,CS,DS,ES,05,FS,#B#1) 


Pour plus de détail on se reportera à l'annexe IT présentant un extrait de LABORDE (1997). 


3.2.6  Evapotranspiration réelle (ETR) 


Cette fonction permet d'évaluer l'évapotranspiration réelle sur un pas de temps quelconque, en connaissant la 
pluie et l'ETP sur ce même pas de temps, la RFU maximale possible et l'état en début de pas de temps de cette 
RFU. On suppose qu'à tout instant le rapport de l'évapotranspiration réelle instantanée à l'évapotranspiration 
potentielle instantanée est égal au rapport de l'état instantané de la RFU à sa capacité maximale. 


RFUmax 100 RFUint RFUinst 
Mois ETP Pluie ETR 15 is ETP Pluie ETR 15 
i 25 45 =HYDROLAB xlalETR($B$1;E2;B3;C3)  =MIN(E2+C3-D3$B$1) 45| 85 515 
2 12 51 =HYDROLAB xlalETR($B$1;E3,B4;,C4) =MIN(E3+C4-D4;$B$1) 51| 8&8. 937 
3 39 75 =HYDROLAB WalETR($B$1,E4,B5;C5) =MIN(E4+C5-D5;$B$1) 75.388 100,0 
4 45 56 =HYDROLABdalETR($B$1;E5;B6;C6) =MIN(ES+C6-D6;$B$1) 56| 45,0. 100,0 
5 86 25 =HYDROLAB.xla!ETR($B$1;E6,B7;C7)  =MIN(E6+C7-D7:$B$1) 25| 65,9. 59,1 
6 123 5 =HYDROLAB xalETR(SBSI:E7.B8.C8) =MIN(E7+C8-D8;$B$1) 5,439 2 
7 156 0 =HYDROLAB xlalETR($B$1;E8;B9;C9) =MIN(E3+C9-D9;$B$1) 0| 159| 42 
8 215 0  =HYDROLAB xalETR(SBS1;E9,B10:C10) =MIN(E9+C10-D10,$B$1) 0, 37 05 
9 165 0  =HYDROLAB dalETR($B$1,E10,B11;C11) =MIN(E10+-C11-D11,$B$1) 0, 04 0.1 
10 57 |15 (=HYDROLAB %alETR(SB#I;E11,B12,C12) =MIN(E11+C12-D12.$B$1) 15| 36 115 
=HYDROLAB-da!ETR($B$1;E12,B13;C13) =MIN(E12+C13-D13;$B$1) 7| 68 837 
=MIN(E13+C14-D14:$B$1 L 152, 178 100.0 


Fig. 83 : Exemple d'utilisation de la fonction ETR 
(à gauche affichage des fonctions, à droite affichage des valeurs numériques) 


Pour plus de détail on se reportera à l'annexe II présentant un extrait de LABORDE (1998). 
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3.3 Fonctions des coordonnées géographiques 


HYDROLAB permet de passer des coordonnées géographiques (longitude et latitude) aux coordonnées 
des projections Lambert fréquemment utilisées. Longitudes et latitudes sont exprimées en ° et dixièmes 
( les longitudes est sont comptées positivement), les coordonnées Lambert sont exprimées en mètres. 

Six systèmes de projection Lambert sont disponibles : France (zone 1 Nord, zone 2 Centre, zone 3 Sud, 
zone 4 Corse) et Algérie (zone Nord Algérie et Sud Algérie). 


3.3.1 Passage des coordonnées géographiques aux Lambert 

Il suffit comme pour les autres fonctions de les faire apparaître dans le menu "Insertion / Fonction / 
Fonctions personnalisées". Les algorithmes nous ont été aimablement fournis par J.P. PIRAT du service de 
géodésie et de nivellement de l'Institut Géographique National, ils sont ceux en vigueur au 15/03/1996. 


Lieux Longitude Latitude | XLambert (m) | Ylambert (m) 

Erwirons de Lille 2,99 50,65, 646172,3115) 328108,5675 

Erwirons de Lyon 4,8 45,7, 791791,5966| 80741,52637 

Erwirons de Marseille 5,2 43,5, 831479,0653| 137367,6587 

Environs d'Ajaccio | 8,8 41,9, 535769,3018| 390860,5389 

Environs d'Alger | 3 37| 526698,523, 410969,6613 

Environs de Ouargla | = 31,9, 764823,6258| 148346,2411 
Lieux Longitude Letitude XLemibert (m) Ylambert (m) 
Environs de Lalle 2,99 50,65 |(=HYDROLAB.xlalXL1($C2,$B2) =HYDROLAB xlalYL1($C2;$B2) 
Environs de Lyon 4,8 457  |=HYDROLABxlalXL2($C3:$B3) =HYDROLAB xlalYL2($C3;$B3) 
Environs de Marseille 5,2 43,5 =HYDROLABxlaIXL3($C4:$B4) =HYDROLAB xlal YL3($C4:$B4) 
Environs d'Ajaccio 8,8 419  =HYDROLABxIaIXLA(SC5:SB5) =HYDROLAB xla! YLA(SC5;$B5) 
Environs d'Alger 3 37 (=HYDROLAB xlalXLNA(SC6,$B6) =HYDROLAB xlal YLNA(SCE,$B6) | 
Environs de Ouægla 15,5 31,9 =HYDROLAB xlalXLSA(SC7 $B7) =HYDROLAB xlal YLSA($C7;$B7) 


Fig. 84 : Exemples d'utilisation des fonctions géographiques 
(en bas affichage des fonctions, en haut affichage des valeurs numériques) 


3.3.2 Passage des coordonnées Lambert aux géographiques 


Le passage inverse des coordonnées Lambert aux coordonnées géographiques, peut être réalisé aisément 
en utilisant le Solveur d'Excel. La figure 85 illustre la démarche qui consiste à minimiser la distance entre deux 
points, le premier connus par ses coordonnées Lambert, le second dépendant de la longitude et de la latitude. 


Point consu Coordormées Géo. | Lambert estimé 
XLambert (m) 646172 Longtude 1 =HYDROLAB xalXL1(D3,D2) 
Ylambest (m) 328108 Latitude 1 Æ=HYDROLAB xial YL1(D3,D2) 
Distance “(B2-E2)2+(B3-E3) 2} OS 


Fig. 85 : Préparation de la feuille 


Fig. 86 : Lancement du Solveur 


Point core Coordannées Géo Lambert estimé 
XLambert (m) 646172 Longitude 2988964,  646099,3267 
Ylambert (m) 328108 Latitude 5064972 328076,2%46 
Distance 79,28839| 


Fig. 87 : Résultats obtenus 
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ANNEXES 


Annexe I : Comblement des lacunes dans des séries pluviométriques 


Bien des études de synthèse sur les hauteurs de pluie annuelles, nécessitent de connaître les précipitations 
sur différentes stations mais sur une même période. C'est le cas par exemple pour l'établissement de cartes 
pluviométriques moyennes interannuelles, où il est vivement recommandé de travailler sur une période 
commune trentenaire. Malheureusement les nv postes pluviométriques n'ont pas tous été observés durant les no 
années de la période de référence et la matrice des données se présente donc généralement ainsi : 


x(1,1) inconnu x(1,j) x(1,nv) 
x(2,1) x(2,2) x(2,j) inconnu 
[x] = x(1,1) x(1,2) X(1,j) x(i,nv) 
inconnu inconnu 
inconnu ds ya 
x(no,1) x(no,2) x(no,j) x(no,nv) 


Si le nombre de poste nv est grand, il est quasiment impossible de "boucher" les trous un par un par des 


techniques de régression. Le choix des variables explicatives est vite inextricable et de plus les erreurs 
accidentelles ou systématiques qui ne manquent pas d'affecter certaines observations, sont ainsi répétées. Nous 
proposons donc une méthode qui permets de combler rapidement et simplement les lacunes avec la partie la plus 


fiable de ce tableau de données. 


Nous calculerons tout d'abord les moyennes expérimentales Mxo(j) et les écart-types expérimentaux 
Sxo(j) des précipitations pour chaque station j et sur les seules années réellement observées : 


[Mxo] = Mxo(1) 
[Sxo] = Sxo(1) 


On peut alors évaluer la matrice des valeurs centrées réduites : uli, j) = 


Mxo(2) 
Sxo(2) 


Mxo(j) 
Sxo(j) 


Mxo(nv) 
Sxo(nv) 


x(i, j)- Mxo(j) 
Sxo(j) 


et calculer pour chaque observation la moyenne des valeurs centrées réduites disponibles : uli) 


u(1,1) inconnu u(1,j) u(1,nv) ui) 
u(2,1)  u(2,2) u(2,j) inconnu u(2) 
[u] = u(i,l)  u(i,2) uj) u(i,nv) ui) 
inconnu inconnu 
mconnu se ERA : 
u(no,1)  u(no,2) u(no,j) u(no,nv) u(no) 


On peut alors dans une première étape remplacer chaque valeur inconnue x(i,j) par une valeur 


Thi) - Mxo(j) 


xegli, j) = So) 


première matrice complète [xo] : xa (i, j) 


x(1,1)  xę&(1,2) x(1,j) x(1,nv) 
x(2,1)  x(2,2) x(2,j) xę&(2,ny) 
[xo] = x(i,1) x(i,2) x(i,j) x(i,nv) 
x@(i+1,1 x@ (i+1,)) 
x@(no-1,2 
x(no,1) x(no,2) X(n0,j) X(n0,nv) 


correspondant à la variable réduite moyenne pour cette observation. On obtient alors une 


35 


Nous effectuons alors une A.C.P. sur cette matrice [xo] et obtenons les projections des variables [az] et 
des observations [c]] sur les k seules premières composantes principales pouvant avoir une signification 
physique : 


c1(1,1) ce1(1,2) … c1(1,k) 
c1(2,1) c1(2,2) … c1(2,k) 
a1(1,1)a1(1,2) .… aj1(1j) … aj(1,nv) e A Le a 
[a] = | a1(2,1) a1(2,2) .… a1(2j) … aj(2,nv) [c11= | c1đ,1) c1(G,2) … cjG,k) 
atly a1(K,2) . ai y ikav 
ci “ae B c1(n00,k) 


Ainsi que les matrices des moyennes et écarts-types : 


[Mx]]= Mx1(1) Mx]j(2) m Mx1() D Mxı (nv) 
[Sx1] = Sx1(1)  Sx1(2) a Sx1(j) ai Sx1(nv) 


Il est alors possible de reconstituer chaque observation manquante x(i,j) par une nouvelle valeur estimée 

xe] (i,j) : 
xe] Gij) = Mx1Q) + Sx1() * { a1(1.j) c1G,1) + a1(2.j) c1G,2) +... + aj(k.j) c1G.k) } 

Cette estimation n'est pas très correcte puisque l'A.C.P. a été effectuée sur une matrice "bouchée" à partir 
de moyennes interannuelles, cependant xe](i,j) est une meilleure estimation que xe, (i, j) puisqu'elle tient compte 
des observations aux autres stations pour cette année j. On peut donc réitérer le processus en remplaçant dans la 
matrice [xo] chaque xe, (i, j) par les xe] (i,j) adaptés. On obtient ainsi une nouvelle matrice [x]] : 


x(1,1)  xe1(1,2) iiy x(1,j) si x(1,nv) 
x(2,1) x(2,2) sis x(2,j) se xe](2,nv) 
[x1] = x(i,1) x(1,2) i X(1,j) yss x(i,nv) 
xe](..,1) = ne xe] (i,j) 
xe] (..,2) 5 a s 
x(no,1) x(no,2) Fe X(n0,j) Le X(no,nv) 


On recommence alors une A.C.P. sur la matrice [x1] permettant d'obtenir de nouvelles matrices [c2] , 
[a2], [Mx2] et [Sx2] d'où l'on tirera de nouvelles estimations xe2(1.j). 


A chaque itération nous 


20 sn : š 

18 modifions les estimations pour 
16 les observations manquantes. A 
14 titre d'exemple, on citera l'étude 
12 —+— Msmumdes écarts de 120 postes pluviométriques 
10 —4— Moyerme des écats| d'Algérie sur une période 
8 —8— ET. des écarts commune de 60 ans où nous 
6 A : 2 
7 avions 1854  stations-années 
2 manquantes sur 7200 (25% de 
0 trous"). Nous avons à chaque 


fois noté la moyenne et l'écart- 
type des modifications, ainsi que 
la plus forte modification en 
valeur absolue. Ces résultats sont 
reportés sur le graphe ci-contre. 


Dif. 0-1 
Dif. 1-2 
Dif. 2-3 
Dif. 3-4 
Dif. 4-5 
Dif. 5-6 
Dif. 6-7 


On y constate que toutes les caractéristiques des écarts d'une itération à l'autre suivent sensiblement une 
exponentielle décroissante. Ici comme dans d'autres tentatives on constate que le processus est stabilisé dès la 
6ème ou 7éme itérations. 
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La difficulté réside dans le choix du nombre de composantes principales à utiliser dans les reconstitutions. 
Nous préconisons bien sûr, de n'utiliser que les composantes ayant une signification physique. La pluie ayant à 
priori une structure spatiale, on peut préjuger que seules les composantes ayant également une structure spatiale 
sont physiquement significatives. Ainsi le choix des composantes à utiliser se basera sur l'analyse du 
variogramme des projections des variables sur les composantes : 


y (h) expérimental y (h) expérimental 
20000 3000 Variogramme Nord-Sud 
15000 = 6000 
Variogramme Nord-Sud Variogramme Est-Ouest 
10000 Variogramme omni-directionnel 4000 Variogramme omni-directionnel 
5000 2000 
Variogramme Est-Ouest 
0 0 
0 50 100 150 200 250 300 0 50 100 150 200 250 300 
Distance en km Distance en km 


Comportement à l'origine de la troisième composante Comportement à l'origine de la quatrième composante 


Dans cet exemple algérien, les trois premières composantes ont une structure spatiale nette (mais 
anisotrope), par contre la quatrième composante est pépitique pure. On ne prendra donc en compte dans les 
reconstitutions que les trois premières composantes. 
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Annexe II : Estimation de l'évapotranspiration potentielle 
par la méthode de Penman et ses dérivées 


IL.1 Estimation de l’E.T.P. par la méthode de Penman originale 


Penman propose d'évaluer l'ETP à partir d'un bilan énergétique simple : 
Rn=A+S+E*L 


Ran représente la radiation nette reçue au sol ; 

A représente le flux de chaleur au bénéfice de l'atmosphère ; 

S le flux de chaleur résultant des échanges thermiques avec le sol ; 
E le flux évaporé 

et L la chaleur latente. 


La formule développée prend la forme suivante : 


ETP ={ £? (1-a) (0.18 + 0,62 À) - -S_ T4(0,56 - 0,08 € ) (0,10 + 0,90 3 } — 
H 59 H = 


5 


¥ ose (ey-e)(1+0,54 V) 


Fr 
Y 


1 + 


ETP  évapotranspiration potentielle en mm/j ; 


Iga radiation solaire directe en l'absence d'atmosphère exprimée en cal/cm?/j | 
a  albédo de la surface évaporante, prise ici pour de la végétation à 0,2 ; 
h durée réelle d'insolation en heures et dixièmes du jour considéré ; 
H durée maximale possible d'insolation en heures et dixièmes pour ce jour ; 
© 
T 


constante de STEFAN-BOLTZMAN soit 1,18 1077 cal/cm2/jour/°K 
température moyenne journalière de l'air sous abri exprimée en degrés Kelvin (T = t + 273 si 
l'on mesure la température t en ° Celsius) ; 
e tension moyenne journalière de la vapeur d'eau mesurée sous abri et exprimée en millibars ; 
F'T pente de la courbe de tension de vapeur saturante pour la température de l'air T ; 


y constante psychrométrique ; 
ew tension maximale possible de la vapeur d'eau, exprimée en millibars, pour la température T ; 


V Vitesse moyenne journalière du vent mesuré à 2 mètres au-dessus de la surface évaporante et 
exprimée en m/s. 
Ces différents paramètres intervenant dans l'évaluation de l'ETP proviennent : 


- de mesures directes sur le terrain pour h, T, et V ; 


- de mesures indirectes pour e, puisque e est estimée à partir de l'humidité relative Hr en %, 
et de la température t en ° Celsius par la relation : 


Es OL HA nait 


- de constantes physiques bien connues 6, F'T, Y et ew : 
o= 1,18 1077 cal/cm2/jour/°K 
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Fr= exp( 19,511 T (es - 3,927) 


(F'T en millibars par ° Celsius, t en ° Celsius, T en ° Kelvin) 


6,149 
T 


y = 0,6605 - 0,826 10-4 z 
(z altitude du lieu en mètres, y en millibar par ° Celsius, 
en fait, on prendra y = Cte =0,66) 


t 
ew = 6,107 10 5 2384 t ) ( ew en millibars, t en ° Celsius) 


- de constantes physiques dépendant de la latitude L du lieu et de la déclinaison D du soleil 
en fonction de la date J . Les formules que nous donnons ici sont issues de R. DURAND 
("Estimation du rayonnement global à partir de la durée d'insolation", Annales 
Agronomiques de l'INRA, 1974). Nous avons comparé les résultats de ces formules avec 
les tables fournies par la F.A.O. (Bulletin N° 24) et y avons constaté des écarts 
généralement inférieurs à 1% ou du même ordre que le nombre de décimales données (5% 
au maximum). 
la déclinaison est donnée par : 

D = 0,0066241 + 0,406149 sin [0,0172029 (J - 81,95) ] 
+ 0,006675 sin [0,0344057 (J - 42,85) ] 
+ 0,003009 sin [0,0516086 (J - 21,42) ] 
+ 0,000149 sin [0,0688115 (J - 17,57) ] 
(D s'exprime en radians et J est le numéro du jour dans l'année, 
de 1 pour le ler janvier à 365 (ou 366) pour le 31 Décembre) 


0,01065 ) 
N cos?(L) - sin? (D) 


( H est en heures, et L représente la latitude en radians ) 


H = 7,6394 { Arc cos[-tg(L) tg(D)] + 


Iga = 914,54 * { sin(L) sin(D) Arc cos[-tg(L) tg(D)] +w cos2(L) - sin?(D) } 
* {1 + 0,033 cos(0,0172 J ) } 


( Iga s'exprime en calories par cm? et par jour) 


` 


- d'une constante physique à évaluer sommairement en fonction de l'état de surface (les 
mesures sur le terrain sont extrêmement rares) : 
nappes d'eau a = 0,05 à 0,07 
cultures vertes a =0,2 
herbes et savane a =0,22 
forêt a =0,11 
roches a =0,16 
sables clairs a =0,26 


IL2 Estimation par la méthode de Penman modifiée F.A.O. 


La formule de Penman a été modifiée par de nombreux auteurs de façon à en simplifier l'écriture ou en 
adapter les paramètres à des régions particulières. Parmi ces formules modifiées on utilise fréquemment au 
Maghreb la formule dite de Penman-FAO que nous donnons ici avec les paramètres recommandés pour les 
régions sèches et arides. 
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ETP = ={ $ (1-a) (0,25 + 0,45 À D --ST(0.56 - 0,08 vE) (0.10 + 0,90 À iyaa 


a (ew-e)(1+ 0,864 V) 


On y retrouve les mêmes variables que citées précédemment, mais seuls quelques paramètres changent : 
(0,25 + 0,45 D au lieu de (0,18 + 0,62 D 
(1 + 0,864 V) au lieu de (1 + 0,54 V) 
IL3 Estimation par la méthode de Penman modifiée par Mc Culloch 
Parmi les modifications apportées à cette formule, citons également celles de Mc Culloch. En effet cette 
modification a été considérée comme bénéfique sous le climat marocain. Cette formule se présente sous la 
forme suivante : 


Fr 
ETP = ={ E (1-a) (0,29 cos(L) + 0,52 d D - £ T0,56 - 0,08 € ) (0,10 + 0,90 D) Ua 


—(ew-e) (1 + 0,864 V) (1 + 


20 30 000) 


On y retrouve les mêmes variables que dans la formule de Penman, mais la latitude L et l'altitude z du 
site (en mètres) sont introduites comme paramètres correcteurs : 


(0,29 cos(L) + 0,52 2) au lieu de (0,18 + 0,62 D 


(1+0,864 V) (1 + 2) au lieu de (1 +0, Fi V) 
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Annexe III : Passage de l'Etp à l'ETR par un bilan tenant compte 
du stress hydrique de la végétation 


La procédure classique de Thornwaite suppose que l'évapotranspiration réelle est égale à 
l'évapotranspiration potentielle, jusqu'à ce que la R.F.U. soit vide. Cependant avant que la R.F.U. soit vide, la 
végétation subi un stress hydrique et diminue son évapotranspiration. Il en résulte que si la part de la R.F.U. 
encore en eau diminue, la végétation n'évaporera pas toute l'E.T.P, mais une part réduite. 

Nous supposerons qu'entre deux instants t] et t2 distants de At il est tombé une quantité de pluie P que 


l'on pourra supposée d'intensité constante : i 


i=— 
b-t 
De même pendant cet intervalle de temps, on supposera que l'évapotranspiration potentielle Ep se 
produit à intensité constante : ep 
sr. 
p = 
t -tı 
Pour passer de l'intensité de l'évapotranspiration potentielle ep à l'intensité de l'évapotranspiration réelle 
er , On supposera que cette dernière est proportionnelle à ep et au rapport de l'état de la R.F.U. r(t), à sa capacité 
maximale R . L'évapotranspiration réelle er(t) est alors variable dans le temps : 


r(t) 
e(t) = DR 
Sur un petit intervalle de temps dt, le bilan en eau de la réserve s'écrit : 
d{r(t)} = i dt - e(t) dt 
Cette relation étant valable tant que r(t)<= R 
On obtient ainsi l'équation différentielle suivante : 


. € 
dfr} =(i Er )dt 
En posant & = ep/R on obtient : 
dé} _ 
r(t) - 1/4 
cette équation s'intègre aisément entre les instants t4 et t : 
[Ln{ r(t) -ia i =[-atl, 


r(t)- 1/0 | _ 
aoa ne en) 


- © dt 


r(t) = i/a + {r(t,) - i/a }e 26t) 


a) cas où r(t2) reste inférieur à R 


r(t2) = i/a + {r(t1) - i/a }e”*(2t1) 


Ep 
r(t2) =E R + {r(ti) E R}e R 


Pour calculer l'évapotranspiration réelle Er entre les instants t1 et t2, il suffit d'intégrer l'équation de e;{t) 


t2 
Er = e(t) dt 
t 
t2 t2 
€p 
Er = R r(t) dt = & r(t) dt 
t ti 


En remplaçant r(t) par son expression calculée plus haut on obtient : 
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t2 
Er= | [i+ {art)- ijen] dt 
t1 
Er= [it] - {rt} [ e26) J2 
Er= i(t-t1) - {r(t) - ia} {e 24241) - 1} 


Ep 
Er=P- {r(tı) -RHE R- 1) 
p 


b) cas où r(t2) deviendrait supérieur R 
On calcule alors le temps t' auquel r(t) atteint R : 
r(t) =i/@ + {r(t,) - je AE) =R 
e2(t-t1) = _R- 1/0 
r(t) - 1/0 
: r(t1) - 1/0 
t= ty + 1/0 Ln tee 


R - i/0 
On calcule alors aisément Er : 


t' t2 
Er -| e(t) af €p dt 
t1 t' 


Er= i(t-t;) - es 4) -1}+ e (tot) 


D La (DEEE - (rh) - iay{e” ds TA - 1}+ ep(t2-t1- 1/0 La =) 
Er= i/a Ln[ OR) . fft) - ia (Ria Bilo- 1}+ eplt-t1)- epla Lopta) 


1/0 


Er = (i/& -e,/0t) pa ia) + ep(t2-t1)- {R - r(t1)} 


c) Mise en œuvre pratique : 
Entre les instants t] et t2 distant de At, il est tombé une hauteur de pluie P, pendant cet intervalle de 


temps At, l'évapotranspiration potentielle est de Ep et à l'instant t], la réserve contenait r(t). 
On calcule tout d'abord : 


Si le résultat est inférieur ou égal à R on aura : 


Ep 
Er=P- (r(t) -RHE R- 1) et r(t2) = à R + {r(t1) À P Rje 
p Ep Ep 


Si le résultat est supérieur à R on aura : 


r(ti) -=> 


r= Ge -RIn(— + Ep- {R - r(t1)} et r(b)=R 


FR 
Ep 
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